Measurement of coherent charge transfer in an adiabatic Cooper pair pump 
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We study adiabatic charge transfer in a superconducting Cooper pair pump, focusing on the 
influence of current measurement on coherence. We investigate the limit where the Josephson 
coupling energy Ej between the various parts of the system is small compared to the Coulomb 
charging energy Ec- In this case the charge transferred in a pumping cycle Qp ~ 2e, the charge 
of one Cooper pair: the main contribution is due to incoherent Cooper pair tunneling. We are 
particularly interested in the quantum correction to Qp, which is due to coherent tunneling of 
pairs across the pump and which depends on the superconducting phase difference ipo between the 
electrodes: 1 — Qp/(2e) ~ (Ej/Ec) cos <po. A measurement of Qp tends to destroy the phase 
coherence. We first study an arbitrary measuring circuit and then specific examples and show that 
coherent Cooper pair transfer can in principle be detected using an inductively shunted ammeter. 



PACS numbers: 74.50.+r, 74.78.Na, 73.23.Hk 



I. INTRODUCTION 

It is well-known that electrons can be transferred 
through a mesoscopic device by means of adiabatic 
changes of system parameters like, e.g., externally ap- 
plied electric or magnetic fields. Since the original work 
of Thouless Q] , several theoretical proposals concerning 
possible physical realizations of this phenomenon, usu- 
ally referred to as parametric pumping, have been put 
forward 0- Depending on the physical mechanism em- 
ployed, these proposals can be divided into two classes. 

In open systems, i.e., systems consisting of several con- 
ducting parts, connected to each other by highly trans- 
missive barriers, parametric pumping can be achieved 
through a periodic modulation of phase of the scatter- 
ing matrix associated to the device [2j- In these propos- 
als, electronic phase coherence plays a fundamental role 
and charge is transferred coherently through the entire 
system. The amount of charge that is transferred per 
period of the modulation is in general a fraction of the 
electronic charge e, which depends on the modulation 
path. In open devices, electron-electron interactions are 
weak, and lead to (small) dephasing corrections to the 
noninteracting result for charge transfer. 

In the opposite limit of closed systems, i.e., sev- 
eral metallic islands connected to each other by ultra- 
small tunnel junctions, parametric pumping of charge 
is achieved by periodic modulation of the so-called 
Coulomb blockade In these proposals, the presence of 
strong Coulomb repulsion between electrons is essential. 
It leads to the quantization of charge on the islands in 
units of e. A periodic modulation of externally applied 
gate- voltages leads to a periodic lifting of Coulomb block- 
ade, which enables the transfer of exactly one electron 



per period through the device. The presence of phase 
coherence in this limit only leads to small corrections to 
the classical result, through coherent higher order charge 
transfer processes known as (in)elastic co-tunneling |fj|. 

Experimental evidence for parametric charge pumping 
in normal metallic systems has been found in both limits 
of open and closed 0, systems. 

The present paper is devoted to parametric pumping 
in a superconducting system Q . It consists of supercon- 
ducting islands, connected to each other by small tunnel 
junctions, and it is operated in the Coulomb blockade 
regime. From this point of view the system is closed, and 
charge transfer is mainly classical, i.e. quantized in units 
of Cooper pair charge 2e per period of the modulation. 
From the other hand phase is defined naturally in a su- 
perconducting system, and due to the Josephson coupling 
between the various parts of the system phase coherence 
tends to be maintained throughout the device. In this 
sense, a superconducting pump behaves as an open sys- 
tem, and the total charge transfer will be characterized 
by significant phase-coherent corrections to the classical 
result |Toj . 

We will be particularly interested in the influence of 
the measurement apparatus of the amount of transferred 
charge on the phase-coherent properties of a supercon- 
ducting charge pump. Measurement on a quantum co- 
herent system is achieved by connecting it to the envi- 
ronment provided by the measuring device. The environ- 
ment in turn can be modeled as a collection of harmonic 
oscillators. Depending on the environment, the quan- 
tum coherence is typically lost in a time which we can 
call the dephasing or the decoherence time once properly 
defined. We focus here on quantum coherent systems 
formed of Josephson tunnel junctions, and their electro- 



2 



magnetic environment. These systems ^2 1X3- , or 
we might already call some of them devices to perform 
quantum logic operations (see review in hold great 

promise in quantum computing because of their potential 
in scalability. 

In a superconducting system the phase difference 
across a junction follows the well-known AC and DC 
Josephson relations. The theory of phase fluctuations 
and environment is well established for electrical circuits 
including small tunnel junctions [l6l UtI ITsj. What we 
do here specifically is the analysis of the back-action of 
the measurement of electrical current on a circuit consist- 
ing of small Josephson junctions. In particular, we show 
that unlike in a standard dissipative ammeter, the phase 
diffusion is limited in a measurement performed by an 
inductively shunted measuring circuit. Measurement of 
tiny currents provides a read-out of Josephson junction 
quantum bits, like the one in the so-called " quantronium'' 
experiment |l4| , or in general when reading out the per- 
sistent current in a flux quantum bit (an RF SQUID loop) 
[l2l ll^ l . As a specific example we focus on a measure- 
ment of current in a double island adiabatic Cooper pair 
pump (CPP), but our conclusions concerning the effect of 
an inductance limited phase diffusion measurement can 
easily be generalized to other Josephson junction circuits 
with straightforward modifications, which account for the 
different topology of the circuit to be measured. 



II. THE MODEL 




CPP 



FIG. f : Double island (a) electron and (b) Cooper pair pump. 
A symmetric pump consists of junctions with capacitance C 
and Josephson coupling energy Ej. Gate capacitance is C g 
for both islands, and the gate voltages are V g \ and V g 2 for 
the first and the second island, respectively. Electron pump 
can be operated with either zero or non-zero bias voltage V 
across; Cooper pair pump has zero bias voltage but it can be 
phase biased by p. In (b) we also show the shorthand symbol 
used in the later figures for a CPP with phase p across. 



the superconducting gap and Eq = 2e 2 /C (C is the junc- 
tion capacitance), the system is in the Coulomb blockade 
regime, and only Cooper pairs are transferred. In this 
regime, the device can be referred to as a CPP, described 
by the Hamiltonian 



H = H C + Hj 



(1) 



An electron pump is a reversible device which provides 
quantized transport of electrical charge upon cyclic op- 
eration of gates connected to it. The simplest variant is 
perhaps a double island pump (see Fig. ^1) in which tiny 
metallic (non-superconducting) or semiconducting grains 
are tunnel coupled to each other and to the surrounding 
electrodes, and capacitively coupled to two gates 0,IH- 
The islands are so small that, due to their tiny capaci- 
tance, equilibrium charge configurations on them are de- 
termined solely by electrostatics on the single electron 
level characterized by the scale of charging energy Ec- 
Applying 7r/2 phase shifted harmonic voltages to the two 
gates at a frequency /, which is so low that it allows the 
system to follow the ground state configuration of charg- 
ing energy at each phase of operation, current through 
the pump equals Jp = ef, based on transport of one elec- 
tron per cycle, on a precision level of a few per cent. A 
more accurate but more complicated device can be built 
by increasing the number of islands and adjacent gates 
in the pump: this approach has been taken to meet with 
requirements in metrology [Tfij . 

If all the electrodes and islands are superconducting 
(see Fig. |IJ>) , charge can be transported not only as sin- 
gle electrons, or rather quasiparticles with single elec- 
tron charge in this case, but also as Cooper pairs |^.ITo|. 
Ideally, in the absence of any bias voltage and at low 
temperatures T such that fc^Tc Ec <C A, where A is 



where 



H c = -E c [(hi - n xl ) 2 + (h 2 - n x2 ) 2 
+ (n 1 - n xl )(h 2 - n x2 )} 



(2) 



is the charging energy of the islands, coupled capacitively 
to each other as well as to two gates 1 and 2. The oper- 
ator ni(2) denotes the number of excess Cooper pairs on 
island 1(2) and n x u X 2) is the gate-charge (in units of 2e) 
which can be tuned by varying the gate voltage V g i^ g2 y 
The Josephson contribution is given by 



Hi = —Ei 



cos(</?/3 + 4>i) 4- cos(i^/3 + 02 - $1) 

+ cos(<^/3 - </> 2 ) • (3) 



Here Ej is the Josephson coupling energy of each junc- 
tion (assumed equal for the three junctions). Through- 
out this paper we will work in the limit Ej <C Ec where 
charging effects dominate. The phase operator 4>\( 2 ) is 
conjugate to the number operator hi^y The phase bias 
<p is the phase difference across the entire pump. As long 
as one can ignore any dynamical environment, (p can be 
considered as a classical variable that can be set to a 
constant offset value tp = (po, for instance, by applying a 
magnetic flux through the circuit of Fig. ^p. 
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FIG. 2: Triangular gating sequence in the n x \, n X 2 plane for 
the pumping cycle as discussed in the text. Hexagons cor- 
respond to regions in which the - schematically indicated - 
classical charge configurations are stable. 



We are interested in the charge transferred through 
the device upon a periodic modulation of the applied 
gate-charges n x i(t) and n X 2{t). Specifically, we assume 
the gate modulation to be adiabatically slow, with a fre- 
quency / <C Ej/hEc- Under this condition the sys- 
tem remains in the ground state throughout the mod- 
ulation if the temperature is low, kpT <C Ej. During 
one period of the modulation, the two-component vec- 
tor n x = (n x i,n x 2) describes a closed path in the corre- 
sponding two-dimensional plane, see, e.g., Fig- EI It can 
be shown that the transferred charge in the ground 
state during one period of the modulation is given by a 
contour integral along the closed path followed by n x , 
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(4) 



Here, we introduced the instantaneous eigenstates |m) 
and energies E m of the corresponding time-dependent 
Hamiltonian QJ, sucn that H(t)\m(t)) = E m (t)\m(t)}. 
The matrix element (I)on corresponds to the current op- 
erator of one of the junctions. For the leftmost junction 
for instance we have 



I = h= Jjsin(<£i + (p/3), 



(5) 



where Ij = 2eEj/H is the Josephson critical current. 
The matrix element is taken between the instantaneous 
ground state |0) and the instantaneous excited state \n). 
An outline of the derivation of Q is given in Appendix A. 

In the limit Ej/Ec — > 0, to leading order tunneling 
events involving more than one junction can be ignored. 
As a result, a periodic gate modulation will lead to a 
transfer of a charge Qp = 2e per cycle through the pump, 
independent of the bias phase ipo, see Appendix lB II This 
is the limit of incoherent Cooper pair tunneling. Ex- 
perimentally, incoherent Cooper pair tunneling is best 
realized by a sufficiently dissipative environment which 
completely randomizes tpo, e.g., by fabricating on-chip 



resistors near the pump, with resistances of the order 
or larger than the resistance quantum for Cooper pairs: 
R K /A = h/{2e) 2 ~ 6.45 kfl Hj. 

Deviations from the classical result are related to 
higher order tunneling processes. In general, pumps in 
the superconducting state are less accurate than their 
normal state counterparts and one fundamental reason is 
the coherence of the superconducting wave function. In 
this paper we are interested in the regime where coherent 
Cooper pair tunneling leads to an important correction 
to the incoherent charge transfer. It provides an inter- 
esting, still unobserved quantum mechanical interference 
correction to the pumped charge (or pumped current in 
continuous cycling) |lCj : 



Q P /2e ~ 1 - 9(Ej/E c ) cos ip . 



(0) 



This result is obtained for a triangular gating, as shown 
in Fig. |21 some details of the derivation are given in 
Appendix IB 21 Although the result is perturbative in 
Ej/Ec, the coherent interference term may be appre- 
ciable. For example for Ej/Eq = 0.1, still rather 
well within the perturbative regime, the variation of the 
pumped current Ip = Qpf, is 1.8 times larger than the 
magnitude of the incoherent current 7i nc = 2e/. The 
coherent contribution to the current can be tuned by ad- 
justing ipo. By assuming a realistic operating frequency 
of / = 100 MHz, 2ef equals 32 pA, and this would pro- 
duce variation in the interference term by 58 pA. De- 
tecting this phase coherent current would not only be 
fundamentally interesting but it would also allow one to 
measure dephasing in a superconducting multi-junction 
circuit coupled to a measuring device. 

We finally note that, in addition to the current Ip = 
Qpf discussed so far, a direct Josephson current 7j os ex- 
ists in response to the phase bias tpo- The Josephson cur- 
rent is maximum, ~ Ij , at the triple point P in Fig. |21 
where the three charge configurations |0,0), |1,0) and 
|0, 1) are degenerate. The gating sequence of interest 
here is away from P: at any given point along the cy- 
cle at most two out of the three aforementioned relevant 
charge configurations are degenerate. At such a degener- 
acy, the Josephson current is smaller, ~ Ij(Ej/Ec). The 
width of the degeneracy is finite, oc Ej/Ec- Hence only 
a fraction of the total cycle contributes to the Joseph- 
son current and 7j os oc Ij(Ej/Ec) 2 siinp . Thus 7j OB can 
be distinguished from Ip , as it depends differently on tpo 
and is independent of the frequency /. 



III. COUPLING TO A MEASUREMENT 
CIRCUIT 

We propose two strategies to measure the pumped 
charge as illustrated in Fig. |31 The first variant, in 
Fig. is a generic inductively coupled measurement 
by a DC SQUID (input coil L, mutual inductance M, 
shunt resistance R) . We call this circuit A in the follow- 
ing. The second circuit, in Fig. Eti represents a directly 
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(b) 




FIG. 3: Two specific examples of inductive measuring cir- 
cuits for current. In (a) the SQUID ammeter is inductively 
coupled to the measured circuit. In (b) current is measured 
by determining the escape probability from the zero voltage 
state of a Josephson junction. The inductance is provided 
by the measuring junction. The phase <p can be adjusted by, 
e.g., applying a flux $ through the loop of the CPP and the 
measuring circuit. 



coupled measurement of current by a Josephson junc- 
tion. In this configuration changes in current of the mea- 
sured circuit are detected as changes of lifetime of a meta- 
stable zero voltage (supercurrent) state of the measuring 
current biased junction |l4l l2l|. This measurement cir- 
cuit will be called B in what follows. The measuring 
junction is biased with a current /; its capacitance Cj 
is the capacitance due to its typically planar geometry. 
When biased near the critical current, I < I c , the junc- 
tion can be characterized by the Josephson inductance 
Lj = v^slrC 1 - I/Ic)^ 1 / 2 '■ This relation conforms with 
the identity between the plasma frequency w p /27r and the 
circuit parameters of the junction: oj p = (LjCj)^ 1 ^ 2 . 

The equivalent circuits corresponding to A and B arc 
shown in Fig. and Fig. Qp, respectively. In (a) we as- 
sume the SQUID to be a pure inductor of inductance L s . 
In both (a) and (b) C v represents the total capacitance 
across the pump circuit, including half of the ground ca- 
pacitance Co of the (symmetric) pump loop. We assume 
that the gate capacitances of the pump are small and 
can be neglected. In (b) the current biased measuring 
junction is replaced by its RCSJ equivalent. 

We do not consider in detail the exact electromagnetic 
environment. But, as will be obvious, the value of the re- 
sistance in shunting the SQUID in A or the junction in B 
does not influence our general conclusion. Furthermore, 
the LC-type lines in the environment, not included in 
the present discussion, further help to decouple the dis- 
sipative environment from the measured circuit. In A 
we can assume that the parallel amplifier impedance and 
in B the dissipative (real) part of the environment are 
included in R. 



FIG. 4: Circuit models of the two measuring devices: in (a) 
we model the SQUID by an inductor with inductance L 3 , in 
(b) the Josephson junction is represented by a parallel LRC 
circuit. 



A few words on the current resolution SI of the two cir- 
cuits are in order. In an optimum measurement by circuit 
A, SI is quantum limited to SI/y/ Af ~ Anh/ L, where 
A/ is the bandwidth of the measurement. This holds 
for perfect coupling to the SQUID. Inserting L = 100 
nH we obtain 6I/y/Af ~ 120 fA and for the bandwidth 
Af = 10 kHz we have SI ~ 12 pA. For faster measure- 
ment one needs higher input inductance: to achieve 10 
pA resolution at 1 MHz, L should exceed 10 fj,H, which 
starts to influence the pump operation itself, as will be 
discussed below. The current resolution of circuit B de- 
pends on optimizing its performance through the measur- 
ing junction characteristics. The ultimate limit of reso- 
lution can be reached when the junction is in the macro- 
scopic quantum tunneling (MQT) limit, where escape to 
a finite voltage state is via tunneling through the wash- 
board potential and thermal escape over the barrier is 
prohibited due to low temperature. Using a trapezoidal 
current pulse of height / and duration At we then have 
the probability P = 1 — exp(— TAt) to escape from the 
zero voltage state, T is the MQT escape rate. Since V de- 
pends on the bias current /, we can use P as a measure 
of the total current through the junction: ideally this 
current is the sum of the bias current and the current to 
be measured, as suggested by Vion et al. 01 > whereby a 
figure of merit of the measurement is SI — SPj (dP/dl) 
giving the smallest resolvable current. Here, SP is the 
resolution in reading P, which for typical averaging is of 
the order of 0.01. By a straightforward analysis we find 
out that there is a trade-off of junction parameters in the 
sense that by reducing SI (by increasing the ratio Cj/I c ) 
we lower at the same rate the cross-over temperature T 
from thermally activated escape to MQT. Using this we 
can find an approximate answer dP/dl ~ {h/2e) / {ksTo). 
This, in turn, with SP = 0.01, and by setting To = 30 
mK, a realistic electron temperature in an experiment, 
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yields 51 ~ 10 pA. The bandwidth of this measurement 
depends on various technical parameters in the setup and 
we will not discuss it here. The comparison of the two 
measurements is limited by the fact that they present 
fundamentally different approaches. Circuit A aims at a 
continuous monitoring of the current at a certain band- 
width. Circuit B in turn tries to grab information of 
the CPP current in a single shot, or by successive cur- 
rent pulses into the escape junction. Both of them can 
achieve the requested current resolution, but the suit- 
ability of the two depends on how they interfere with the 
measured CPP and whether one can initialize the CPP 
controllably. One of the advantages of circuit B is that 
the phase fluctuations are smaller owing to the smaller 
inductance of this detector. 

In the presence of the measuring circuit the bias phase 
difference tp acquires its own dynamics. As a result, the 
phase will contain a fluctuating part, ip(t) = (po + 4>(t), 
and this will in general modify the pumped charge. Be- 
fore studying this modification in detail for the set-ups 
of Fig. |3 we analyze the time-dependence of phase fluc- 
tuations induced by the measuring circuits. Fluctuation- 
dissipation theorem implies that the variance of the fluc- 
tuations obeys the relation 



! (i)>=8 



duj 5RcZ(cj) 
Rk 



coth( )[1 — cos(wt)l. 



(7) 

In quantifying the phase fluctuations, we thus need to 
know the real (dissipative) part of the impedance seen by 
the measured circuit, $teZ(w), where w is the (angular) 
frequency. 

Circuit A — In this case we obtain 



$teZ A (t 



M 2 Ruj 2 



R 2 {1 - LC p to 2 ) 2 + lv 2 L 2 [1 - (1 - 7 )LC> 2 



(8) 



with 7 = M 2 1 (L S L). For typical experimental values this 
represents a narrow resonance peak at to ~ lolCi where 
lulc = (LCp) -1 ' 2 is the natural LC resonance (angular) 
frequency of the CPP circuit alone. We first study the 
limit 7 <C 1, i.e., weak coupling of the measuring circuit. 
The circuit can then behave in two different ways de- 
pending on whether the parameter a 2 = (to lcL s 1 1 R) 2 
is either « 1 or » 1. In the first limit the reso- 
nant frequency is ojlc and the width of the peak is 
5lo/cj L c =s {M/L) 2 JlJC~ p /(2R) = 7 a/2 < 1. The 
product of the width and height of the resonance peak 
is $leZ(Lo) max 5u ~ tOLCy/L/Cp, independent of R. In 
the second limit, a 2 3> 1, which corresponds to most 
experimental conditions, the measuring circuit decreases 
the effective inductance of the resonant circuit and the 
resonant frequency (maximum of $teZ A (to)) is shifted 
very slightly upwards to tOLc/y/l — 7- The width of 
the peak is now Slo/uj L c — {M/L S ) 2 R/ (2y/L/C p ) = 
7/ (2a) -c 1. The product of the width and height is 
again $teZ(uj) rnax 8uj ~ UJLC\/L/Cp, independent of R. 



From Eqs. Q and JSJ we see that in both limits the 
phase fluctuations do not diverge with time. Instead 
we have initial oscillations which dephase and asymp- 
totically the fluctuations tend to a constant expecta- 
tion value. This behavior is not so unexpected since 
the measured circuit forms essentially a SQUID loop, 
where the average voltage over the junctions vanishes. 
The decay time of coherent phase oscillations turns out 
to be ^ec, A = LC P R/(^L S ) when a 2 < 1 and t dec , A = 
L s /( 7 i?) in the opposite limit. We obtain the following 
expressions for the variance in the short and long time 
limits, respectively, 



4.(^ 7 ^)coth(^ C 



Rk 

[1 - COs(u! LC t)} 

1 



2k B T' 

if t <C t de c,A 
if t » t deCt A 



(9) 



In the opposite limit of strong coupling of the measuring 
circuit, 7 — » 1, we essentially approach circuit B to be 
discussed below. 

Circuit B — In this case we have 



sReZ B (Lu) = 



L 2 Ruo 2 



R 2 (l - LjC*,uj 



2\2 



,2 T 2 



(10) 



with C} = Cj 



C p . In this case the low damping 



limit corresponds to Suj/lulc — y/Lj/Cj/R -C 1, where 
ujlc = (-kjCj) 2 is again the LC resonance (angu- 
lar) frequency of the CPP circuit now with the mea- 
suring junction. In this limit the maximum of the real 
impedance is ^ReZB(to) ma x — R and the product of height 
and width is ^ReZB{to) max Su> ~ lulc\/Lj/Cj, again in- 
dependent of R. 

The coherence of oscillations in the underdamped case, 
Slu/ujlc *C 1 dies in a time tdec,B = 2RC* ; , and as a 
result, the asymptotics of the variance for short and long 
times are given by 



(A*)>^47T( 



Rk 

[1 - cos(w LC i)] 
1 



if t < t d 

ec.B 

if * > t d ec.B 



(ii) 



It is interesting to see how the resonance behavior char- 
acteristic for circuits A and B transforms into an essen- 
tially diffusive one when the quality factor of the reso- 
nance in the dissipative LC circuit gets lower. Let us 
consider circuit B as the example. In the overdamped 
case, Suj/ujlc S> 1, the phase fluctuations are diffusive 
when t -C td ec - In particular, when k B T Hujlc, we 
have 



{(j> 2 (t)) ~ 8irk B TRt/(hR K ), 



(12) 



like in a purely dissipative environment 22]. At t ~ 
8lo/w\ c the phase fluctuations level off at the value given 
by the long-time, high-temperature limit in Eq. (|ll|l . 



{cj> 2 {t)) ~8nk B TLj/(hR K ). 



(13) 
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FIG. 5: (a) Numerically calculated ((f> ) for circuit B. In (b) 
we show $teZ(ui). The parameter values for these graphs 
are: Lj = 0.1 nH, C} = 0.1 pF, R = 1 kfi, T = 50 mK, 
whereby 5lu/lo lc = 0.03162 and lu lc = 3.162-10 11 s" 1 . (c) 
and (d): demonstration of the influence of the overdamped 
inductive environment. For both graphs we used the follow- 
ing parameters: Slu/ulc = 10, hu>LC /ksT = 0.1, R — 5 Q. 
The predictions of Eqs. 11121 and I13H are shown by the dot- 
ted and the dashed lines, respectively, and they intersect at 
ULct = Slu/ulc- 



In the limit Suj/lolc 3> 1 the circuit attains features on 
one hand of an inductance limited phase diffusion mea- 
surement and on the other hand of a dissipatively de- 
phasing environment. 

Illustrations of the results of this section for model cir- 
cuit B are shown in Fig. [5] In Fig. [5^ we see the results 
of numerically calculated (4> 2 (t)) for typical parameter 
values, and Fig. 03 shows $ieZ(u>) with the same param- 
eter values. Fig. [3]; demonstrates the crossover suggested 
by Eqs. (|12l) and (| 1 31) in the dissipative case; the corre- 
sponding $teZ(w) is shown in Fig. [SJi. 

The analysis of the variance of phase fluctuations per- 
formed in this section indicates that the measuring cir- 
cuits A and B, unlike a purely resistive environment, 
will not lead to a complete suppression of coherence of 
the CPP: the phase diffusion is essentially limited by 
the inductance. In the next section, we will analyze 
the influence of these fluctuations on the amount of (co- 
herently) transferred charge through the pump. As we 
will sec, fluctuations renormalize not only parameters of 
the pump, but may also induce (coherent) higher order 
charge transfer, like co-tunneling. Both effects lead to 
additional corrections to the pumped charge. 



IV. INFLUENCE OF THE MEASURING 
CIRCUIT ON THE TRANSFERRED CHARGE 

A convenient way to discuss the influence of the mea- 
suring apparatus on the charge transferred during one 
pumping cycle is to formulate the problem in the frame- 



work of the so-called effective action approach. This ap- 
proach enables one to obtain the partition function of 
the pump together with its measuring environment as a 
path integral. For the case of an inductive environment, 
the degrees of freedom of the environment can be inte- 
grated out approximately, and we obtain an effective ac- 
tion that describes the low-energy behavior of the pump, 
in the presence of the measuring circuit. From this ef- 
fective action, we will obtain an effective Hamiltonian 
for the pump, which is essentially the Hamiltonian Q 
up to two modifications: (i) the Josephson energy Ej 
in 121 should be replaced by a renormalized value Ej, c g 
and (ii) a term H{ n d should be added, which accounts 
for correlated tunneling events induced by the measuring 
circuit. Both modifications will lead to corrections to the 
result ©. 



A. Effective action 

We start our analysis by considering the total equilib- 
rium partition function of the pump together with the 
measuring circuit, 



Z,, 



where the Euclidean action S of the system is given by 
S = Sc + Sj + Sbath- Here, Sc and Sj are the actions 
associated to the Hamiltonians He and Hj defined in 
Eqs. (0) and ©, respectively. Specifically, the action Sj 
is given by 



E 



3 



dr cos[<50(T)j + 0(t)/3], 



where we introduced the phase differences across each 
junction in the CPP: S<f>i = <^o/3 + 4>\, 5(f>2 = <po/3 + 
02-0i, and S(f> 3 = tpo/3 - 4>%, The third term, Sbath, is 
the action describing the measuring circuit. Due to this 
contribution, the phase bias ip has acquired a dynamical 
part, 4>{t), such that f{r) — tpy + 4>{t). The action Sbath 
can be written as 



5, 



bath 



dTdT'^D^T-T')^'), (14) 







where (3 = 1/T (in this subsection we use units such 
that h — ks — 1). The kernel D$(t) is related to 
the impedance of the measuring circuit [2j|, D^t) = 
TJ2n D d^n)e- iUnT , where 



1 



4e 2 Z(i\u n \) 



(15) 



Here, w n is a bosonic Matsubara frequency, u> n — 2nnT. 
It is possible to get an effective action depending on the 
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phases <fix (r) and </> 2 (r) of the islands only. We write the where we introduced the renormalized Josephson cou- 
total partition function as pling 

Z tnt ~ I VfaVfoe- 8 **. Ej, cS = Eje~ G ^ T = ^ 18 , (17) 

with the phase-phase correlation function 

G <A (r)=T^D- 1 (o; n )e--« T . (18) 



To second order in Ej, S e g is found by a re- 
exponentiation of the averages (Sj) and (Sj) over the 
action Sbath- For the average of Sj we find 





(Sj) = E JiCS ^T J dr cos S&ir), (16) For the average of Sj we obtain 



P 2 3 



(Sj)^-T^- I rfnrfT 2 {cos^ i (ri)cos^(r 2 ) £ ^(n-^/g + sin( 5^( ri ) sin^(r 2 ) £ CTe ^,(r 1 -T,)/9 } _ 

- 1 ^ LI (T = ±l 

! [191 



a=±l 



We see that at this order an interaction term appears, 
whose kernel depends on the correlation function G$(t). 
As already discussed in the previous section, fluctuations 
in the bias phase are bound for the particular measuring 
environments we are considering. In the limit of small 
phase fluctuations, the correlations ~ (<^(r)0) remain 
small, allowing us to expand the exponentials in (|19f) with 
respect to G^. Therefore we can write 



{Sj)=E J cff ^r J d,TldT2[cOS 5(t>i(Ti) COS 5(f) j(T 2 ) 

•sin (J^i(ri)G^(ri - r 2 ) sm<%(r 2 )]. (20) 



*,3'=1 b 
1 . „ 

~9 



We are now in a position to obtain the effective action 

SeS, 



where 



S e fi = Sc + Sj e s + Sind j 



3 r 

Sj,eS = E JtCS *^2 / dr cos 6(f>i(r) 



(21) 



and 



E f 
Sind = ^— / dridr 2 x 

sin<5</) i (Ti)G'0(Ti - r 2 ) sin^j(r 2 ). (22) 

The influence of the measuring circuit is indeed two-fold: 
(i) In the effective action, the bare Josephson coupling 
energy Ej has been replaced by an effective, renormal- 
ized value, Ej tC ff; (ii) The measuring circuit correlates 



tunneling events at different junctions, as is seen from 
Eq. 1221) . Physically, a tunneling event occurring at time 
Ti in junction i virtually excites the measuring environ- 
ment, which, upon relaxation, causes a second tunneling 
event in junction j at time r 2 . In other words, the en- 
vironment induces higher order tunneling events in the 
pump of the co-tunneling type. Clearly, both (i) and (ii) 
will modify the charge transferred through the pump, as 
we will discuss in more detail below. 



B. Renormalized Josephson coupling 

In order to evaluate the renormalized Josephson cou- 
pling energy E j je g as given by Eq. I|17l) , we need to calcu- 
late the correlation function G$, Eq. (|18fl. at time r = 0. 
This implies calculating the sum over all Matsubara fre- 
quencies of the quantity Z(i|w„|)/|a; n |, which involves the 
total, imaginary-time response of the circuit, Z(z|w„|). 
For the circuits A and B of interest here we find 



Z A (i\uj\ 
RL + (1 — j)L s L\oj\ 



R(l + LC p \uj\ 2 ) 



and 



Z B (i\u>\)/\u\ 



L S M[1 + (1- 7 )LC P M S 



RLj 

fl(l + LjC}H 2 ) + HLj' 



(23) 



(24) 



respectively. Note that, upon analytic continuation 
iu> n — ► uj + iO, the real part of Z A ,B(i\^>n\) coincides with 
the results JHJ and (|10|l . respectively. 

From a closer inspection of Eqs. i|23|) and l|24|) we see 
that we can obtain Zb from Za upon taking the limit 
7 — > 1 in (|23|) . thereby replacing L and L s by Lj. We 
will, therefore, focus on circuit A in what follows; the 
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1 . 



FIG. 6: Renormalization parameter 770 for circuit B in units 
ojlcL/ Rk as a function of a and 7. 

results obtained can be used to analyze circuit B after 
taking proper limits. 

In the zero temperature limit T — > 0, the sum over 
io n in Eq. 1|18|) can be replaced by an integral. In this 
limit, the renormalized Josephson coupling energy can 
be written as -Ej, c ff = Eje~ no , where 

00 

_ 4uj lc L r 1 + a(l - -f)x 

V °~ QRk J X 1 + x 2 + ax[l + (1 - 7 )a; 2 ]' 


We introduced the dimensionless integration variable 
x = oj/ujlc and used the parameter a = u!lcL s /R intro- 
duced in Section ITTT1 Note that the integral is always fi- 
nite: the integrand is well-behaved, both for small values 
of x, reflecting the fact that low- frequency phase fluctua- 
tions are limited by the presence of an inductance in the 
circuit, and for large values of x, reflecting the natural 
high frequency cut-off of the fluctuations provided by the 
presence of capacitances. In Fig. [H] we plotted 770 as a 
function of the dimensionless quantities a and 7. Let us 
use some realistic parameter values to estimate the effect. 
For circuit A, a ~ 1 and 7 -C 1 ; hence we have 770 < 0.1 
and Ej te s > 0.9Ej. For circuit B, a ~ 1CT 2 and 7=1; 
hence we have rjo < 10~ 3 and Ej iG g ~ Ej. We conclude 
that the renormalization of Ej is always weak. This re- 
sult confirms the analysis of Section ITTT1 for small values 
of a or 7 the parameter 770 indeed equals the asymptotic 
values for (</> 2 )/18, see Eqs. © and l(TTjl. Thus, for suit- 
ably chosen parameters, the measuring circuit A (as well 
as B) will not suppress the coherent coupling between 
various parts of the pump. 

The results presented so far are valid at zero temper- 
ature, T = 0, and in the absence of a gate modulation, 
/ = 0. Finite temperature introduces a correction to 
770 due to the summation over Matsubara frequencies in 
Eq. 1|17|) . As long as the inductance L of the pump is 
small compared to the inductance H/ (2elj) of the Joseph- 
son junctions, this correction can be neglected at temper- 
atures of interest here kpT <C Ej. The effect of a finite 
frequency / can also be ignored at this point, since we 
are working in the leading order in the adiabatic approx- 
imation. 



The above results, obtained for an inductive environ- 
ment, should be contrasted with the effect of a purely 
resistive measuring circuit, characterized by a resistance 
Rq. In this case, the low- frequency impedance is con- 
stant, Z ~ Rq, and the sum over Matsubara frequen- 
cies in Eq. (|18fl diverges logarithmically. This reflects 
the fact that the low-frequency fluctuations of the phase 
are unbound. As a result, the Josephson coupling will 
be significantly renormalized. Indeed, in the case of a 
very resistive environment, characterized by a large di- 
mensionless resistance ratio AR^/9Rk 3> 1, the renor- 
malized Josephson coupling at T — is given by £?j >e ff ~ 
EjifRoC p ) 4R °/ 9RK , which depends explicitly on the op- 
crating frequency / of the pump. This means that the 
Josephson energy is completely suppressed in the absence 
of a gate modulation, / = 0. For finite values of /, the 
pump will be characterized by a finite Josephson coupling 
energy; however, its value will be small, Ej eS /Ec <C hf, 
for any reasonable operating frequency. We thus con- 
clude that it is practically impossible to operate the 
pump in the adiabatic regime if Rq is large. In the oppo- 
site limit of a weakly resistive environment, ARq/QRk <S 
1, one finds Ej te $ ~ Ej(EjR C p /h) 4R ^ 9RK < Ej: the 
renormalization of the Josephson coupling is limited. We 
conclude that, at least in principle, the device can be op- 
erated coherently in the adiabatic regime in the presence 
of weak dissipation, as long as Rq is sufficiently small. 

C. Calculation of the transferred charge 

We finally turn to the actual calculation of the charge 
transferred per cycle. In order to use Eq. which 
involves instantaneous eigenstates of a Hamiltonian, it 
is convenient to transform the effective action S e g into 
an effective Hamiltonian, H c r = He + Hj c g + Hind- 
With the help of Eq. \2l\i we see that Hj c g can be ob- 
tained directly from ©, replacing Ej by the renormal- 
ized value Ej <e f[. The Hamiltonian -ffi n d is obtained from 
Si n d, which depends on the correlation function G$. The 
fact that is in general non-local in time complicates 
matters. However, for the measuring circuit A, the dy- 
namics at low frequencies is dominated by the inductance 
L (for circuit B it is Lj). Provided that L is small, such 
that L <C H/i2eIj), the relevant correlations arc local 
in time and G^t) ~ (8ttL/Rk)S(t). In other words, 
the two correlated tunneling events induced by the mea- 
suring environment are instantaneous on the slow time 
scale characterizing the junction dynamics. Substituting 
the above local form for the correlator into (1221) and 
performing one of the integrations over imaginary time 
we see that the resulting Si n d can be transformed into 

3 

ind = 18~ (h/2e) 2 ^ Sm ^ Sm ^ J ' ( ' 

The calculation of Qp with 10} with the Hamiltonian 
-ffeff for a triangular gating sequence is straightforward. 



9 



Ignoring Hind, one obtains the result ©, with Ej re- 
placed by its renormalized value Ej tC g. If -ff; n d is taken 
into account, a correction is found due to coherent higher 
order tunneling events induced by the inductive environ- 
ment. Referring the reader to Appendix IB 31 for details, 
we only present the final result: 

Qp^e^l-^^-.-^-^jcos^ . (26) 

Taking Ej/Ec = 0.1 with Eq ~ IK, the renormalization 
of Ej is negligible. With L = lOOnH (circuit A), the 
second correction Ej^gL/{6(h/2e) 2 ] ~ 0.2, whereas for 
circuit B this would be Ej, cff L/{6(h/2e) 2 } ~ 2 • 10~ 4 . 



V. DISCUSSION 

The CPP has a few properties that make it an attrac- 
tive object of investigation both theoretically and exper- 
imentally. The mechanism for deviations of the pumped 
charge from the quantized value are interesting in their 
own right. These include Landau-Zener (LZ) band cross- 
ing, quasiparticle current, and last but not least the topic 
of the current manuscript, the coherent quantum inter- 
ference in Cooper pair transport. They all can be avoided 
at least in principle by using a favorable operation mode: 
LZ crossing can be minimized by operating the pump 
at low enough frequency, / <C Ej/(HEc), quasiparticle 
current by low temperature and careful filtering of the 
sample, and coherence induced corrections by dissipative 
environment of the CPP. 

Quantized charge transport is being investigated 
largely because it may eventually fulfill the requirements 
of providing a modern standard for electrical current. At 
the same time it would close the metrological triangle of 
electrical quantities: by now both the voltage and resis- 
tance have their quantum standards, Josephson voltage 
and quantum Hall resistance, respectively, and by deter- 
mining the current using I — qf would test the validity 
of Ohm's law with the same value of Planck's constant in 
the two. But there are serious problems to be solved be- 
fore I can be determined at sufficient absolute accuracy, 
error rate being smaller than 10 -7 and output current 
on the nA scale. The foremost problems are either the 
very small current obtained in the tunnel junction based 
pumps, or the errors in the number of electrons carried in 
the moving quantum dots in the semiconductor pumps. 

Cooper pair pump may solve these problems eventu- 
ally. The current limitation can perhaps be lifted by 
using higher values of Ej employing Josephson junctions 
made of a superconductor with higher T c than that of 
aluminum (1 K). An obvious candidate is niobium, where 
T c ~ 9 K, whereby yielding an almost one order of mag- 
nitude enhancement in Ej cx T c . Fabricating small junc- 
tions using Nb has, however, turned out to be a challenge, 
which has not been fully solved yet [24l |25| . The second 
problem is how to suppress errors due to the quantum 



interference in the CPP without disturbing the opera- 
tion of the pump otherwise. As an example, by insert- 
ing a highly dissipative termination to the pump, heat- 
ing becomes a problem especially at the desired higher 
throughput currents. Therefore understanding the na- 
ture of this interference is of importance in optimizing 
the operation of the CPP. One alternative to dissipative 
environment is to employ pumps with larger number of 
junctions, N, since the interference correction is propor- 
tional to (Ej/E C ) N - 2 . 
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APPENDIX A: DERIVATION OF EQ. © 

The Hamiltonian H, Eq. QJ, depends explicitly on the 
gate-charges n x i and n X 2 through the charging term He, 
Eq. J5J. As a result of the periodic modulation of the 
vector ri x — (n x \, n X 2) as a function of time, the Hamil- 
tonian becomes time-dependent. If the modulation is adi- 
abatically slow, Schrodinger's equation ihd\ip m (t)) / dt = 
H (t)\tp m (t)} is satisfied by an adiabatic state of the form 

|V> m (f)> = e-^'^^/V^Mf)), 

where \m(t)} is an instantaneous eigenstate of Hit) with 
energy E m (t) and j m is Berry's phase p6|. 

The expectation value of any time-independent hermi- 
tian operator O in the state \ip m (t)} is in general time- 
dependent and given by (0(t)) = (i/) m (t)\0\i/) m (t)) = 
(m(t)\0\m(i)) . Specifically, we see that the dynamics 
of (0(t)) is governed by the instantaneous eigenstates. 
Using the fact that \rh) = n x ■ <9gjra), where a = da/dt, 
we obtain 

(6) = 23?e (m\6\n)(n\dfi x m) ■ k x . (Al) 

Note that the term n — m can be excluded from the sum 
in (|A1|) . as (n\dff, x n) is purely imaginary. This can be 
seen using the normalization condition (n|n) = 1 from 
which we obtain dn x (n\n) = = 2!fte (n\df{ x n). 

We are interested in the total charge Qp transferred 
through the pump during one cycle. Since the mod- 
ulation is periodic, Qp is the total charge transferred 
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through any of the junctions during the cycle, e.g. the 
leftmost one, Q\ . Let us consider the matrix elements of 
the current passing through the left junction, (m\Ii\n), 
see Eq. |JS}, 

(m\Ii\n) = Um\[H,Qi]\n) = UE m - E n )(m\Qi\n). 

(A2) 

Applying the result IjAlfl to the charge operator Qi and 
using HA2|1 . we obtain 



(Ql) = 2ftQta 



(m\Ii\n) 



(n| dn x m) ■ ri x . (A3) 



Therefore, the total charge passing through the pump in 
the state \ip m } is given by 



Qp= / dt{Qi 

J cycle 



2fr3m 



E 



(h)r. 



(nld^m) ■ dn x 



.(A4) 



which is the result Eq. (0J upon setting \m) = |0). 



APPENDIX B: TRANSFERRED CHARGE FOR 
TRIANGULAR GATING 

In this Appendix we outline the perturbative calcula- 
tion of Qp with Eq. (@J , for a triangular gating sequence 
as in Fig. This sequence corresponds to a contour con- 
sisting of three linear segments, (1), (2), and (3), in which 
respectively the gate-charge vector n x is changed adi- 
abatically from (0,0) to (1,0), increasing n x ±; then from 
(1,0) to (0,1), simultaneously decreasing n x \ and increas- 
ing n X 2', finally from (1,0) to (0,0), decreasing n X 2- In this 
Appendix, when evaluating Qp, we will present explicit 
calculations for segment (1) from (0,0) to (1,0) only, i.e. 
we will calculate 



Q 



(i) 



2 Krai 



E 



(Il)on 
Eq — E T , 



(n\d nxl 0)dn x i 



(Bl) 



Results for the segments (2) and (3) will be simply stated, 
their calculation being essentially analogous to the one 
for segment (1). 

Since Ej <C Eq , only a limited number of charge states 
need to be taken into account in the calculation. This en- 
ables one to proceed perturbatively; below we will con- 
sider various contributions to Qp that arise in different 
orders of perturbation theory. 



1. Incoherent contribution 

In leading order, the only relevant charge states for 
the segment (1) of interest here are the eigenstates |0,0) 



and |1, 0) of He, Eq. J5J, which are mixed coherently by 
Hj, Eq. 10. This is analogous to the coherent mixing of 
charge states in a single Cooper pair box, see 0. The 
two lowest instantaneous eigenstates are therefore 



|0) =a|0,0)+e iVo/3 6|l,0), 
|1) =6|0,0)-e^ o/3 a|l,0), 



(B2) 
(B3) 



with an energy difference E% — Eq = Ej Vl + e 2 ■ Here 
a 2 = 1 - b 2 = (1/2)(1 + e/yT+lF) such that ab = 
l/(2Vl + e 2 )- The dependence on the parameter n x \ en- 
ters through e = (2^/3^) (1 — 2^1). Higher states can 
be ignored as they are separated in energy by an amount 

- E. ■ 

We proceed by evaluating the various terms appearing 
in IjBlfl . In order to obtain the matrix elements (h)oi for 
current through the leftmost junction, it is convenient to 
use the charge representation of JSJ. Putting (p = ipo we 
find 

* = % E (e itpo/3 \n 1 ,n 2 )(n 1 + l,n 2 \ 



"1 ,112 



^Ini + l.TteXm,^!). (B4) 



A direct calculation, using the above decomposition of |0) 
and |1) into charge states, then yields (//)oi = 
Similarly, the decomposition can be used to calculate 
(l|<9 nxl 0) = bd nxl a — ad ri:cl b. Finally, using the equal- 
ity Ei — Eq = Ej/2ab, we obtain, in leading order, 

1 

Q$ ~ 4e J ab(ad nxl b - bd n:cl a)dri x i 



Be / a(l ~ a )da = 2e. 





(B5) 



In this order, the contributions from segments (2) and 
(3) are zero, and we conclude that Qp — 2e. Indeed, in 
the absence of higher order tunneling processes, the cur- 
rent through the leftmost junction is not affected by the 
coherent mixing of charge states |1,0) and |0, 1) along 
segment (2) or of |0, 1) and |0, 0) along segment (3). Be- 
low, we will include higher order tunneling which leads 
to corrections to these results of order Ej/Ec- 



2. Coherent correction 

We now take into account the charge states that Hj 
mixes into |0) and |1) by second order tunneling pro- 
cesses. For the segment (1) of interest, these are the 
charge states |0,1) and |1,— 1). Straightforward second 
order perturbation theory yields the correction to the 
states |0) and |1), 



<5|0> = c\l) and 5\l) = c*|0), 



(B6) 
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where 

c = (3Ej/E c )ab(b 2 e lipo - a 2 e - lipo ). (B7) 

The energies Eq and E\ are also renormalized, yielding 
a correction to the energy difference 

8{E X - E Q ) = {6Ej/E c )ab cos ip . (B8) 

Note the appearance in this order of a dependence of the 
corrections on the bias phase ip . Along segment (1), 
phase coherence through the entire pump is established 
by coherent mixing at the leftmost junction in combina- 
tion with second order tunneling through the middle and 
rightmost junctions. 

The matrix elements of Ii are not affected by the cor- 
rections. Using HB6|) we obtain S(Ii)oi — (=*[(/; )oo + 
(iz)ii], which vanishes since the instantaneous eigen- 
states carry no current: (ij)oo = (Il)u — 0, as can 
be seen easily from 1)B2[L (|B3|I and l|B4|) . However, 
(5(l|9„ xl 0) is non- vanishing, <5(l|<9 raxl O) = <9„ xl c. The re- 
sulting correction to can be written as the integral 

n P = 7^-cos<^o / dn x i{24(a6) 3 (6<9„ a - a<9„ xl 6) 
2e E c J 
o 

-6ab[{b 3 - 3a 2 6)<9„ xl a- (a 3 - 3 2 a)d nxl b]}, (B9) 

where the first term stems from the energy renormal- 
ization (|B8|I and the second term from the correction 
<5(l|<9„ xl O), calculated taking a derivative in i|B7|) . Direct 
integration yields 

SQp-'/2e — — 3^- cosipo. 
Ec 

A similar calculation for segments (2) and (3) yields the 
same result. In other words, the total correction to Qp is 
given by SQp = —2e(9Ej/ Ec) cos^Oi m agreement with 
Eq. 0. 



3. Correction due to inductive coupling 

We finally consider the corrections to Qp associated 
with the Hamiltonian Hi u d, Eq. Q25J1 ■ which depends 
quadratically on Ej. Hence, first order perturbation the- 
ory in -ffind yields a correction to Qp which is linear in 
Ej. We first find the corrections to the lowest instanta- 
neous eigenstates 

<5|0) =d|l) and£|l) = -d*\l), (BIO) 

where 

d = {Ej. cS L/[18(h/2e) 2 }}ab(a 2 e~ iV0 - b 2 e lv "). (Bll) 

The energy difference is renormalized as well, 

5{E 1 - E ) = -{E 2 cff L/[9(h/2e) 2 }}abcos VQ . (B12) 
Comparing these results with the corresponding ones 
found in|Bl Eqs. fEU), ifBTjl . and |B8|l . we conclude that 
the subsequent calculation will be completely equivalent 
to the one performed in IB 21 Indeed, the matrix elements 
of I\ are not affected by the corrections, and the contri- 
bution from i5(l|9„ xl 0) is equal to d Uxl d. As a result, the 
correction <5Qp| nd can be presented as an integral similar 
to the one in (|B9|) . Segments (2) and (3) give the same 
contribution and the total correction to Qp due to Hi n d 
is therefore given by 

SQp M d = 2e{E, LcS L/[6(h/2e) 2 }}cosp , 

in agreement with Eq. I|26|l . The result is perturbative, 
the parameter Ej^ c gL/(h/2e) 2 must be small compared 
to unity. This is in agreement with the condition L <^ 
(h/2elj) stated in the main text for the derivation of the 
Hamiltonian Hind- 
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